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ABSTRACT 

In this paper the generalisation of previous author’s formulation of fractional continuum mechanics to 
the case of anisotropic non-locality is presented. The considerations include the review of competitive 
formulations available in literature. The overall concept bases on the fractional deformation gradient 
which is non-local, as a consequence of fractional derivative definition. The main advantage of the pro¬ 
posed formulation is its analogical structure to the general framework of classical continuum mechanics. 
In this sense, it allows, to give similar physical and geometrical meaning of introduced objects. 


1 Introduction 

Reliable modelling of heterogeneous materials, ranging from marco- to micro- and nanoscale of obser¬ 
vation, in terms of continuum mechanics concept, needs non-local formulations [1]. The first articles in 
this area were released in the 1960s, and dates back to the scientists, such as Toupin [2], Mindlin [3], 
Eringen [4], Dillon [5], Dafalias [6], Bazant [7], Maugin [8], Aifantis [9], Fleck and Hutchinson [10], 
Regardless of the details of specific formulation, it is common that non-local model introduces charac¬ 
teristic length (or time) which is inherent to the inner material structure. This new material parameter 
can be directly measured from experiment or can be determined from experiments via inverse analysis 
or derived theoretically from micromechanics [11, 12, 13, 14, 15, 16], 

Nowadays there are many concepts dealing with non-local formulations, like general non-local theories 
[17, 18], strain-gradient theories [2, 19], micropolar theories [20, 21] or theories of material surfaces 
[22]. Nevertheless, due to development of new materials, and constant miniaturization of e.g. electronic 
or medical devices, the continued progress in this subject is desirable. In some sense, as a response to 
this demand the new branch of investigations was initiated in early 2000’s by Klimek [23], Vazquez 
[24], Lazopoulos [25], namely the one which utilises fractional calculus. 

An important class on non-local models are those dealing with anisotropic non-locality - especially 
the one combined with damage description. Such formulations introduce the possibility of analysis of 
materials showing scale effect but depending on direction. Let us mention herein, the papers in this 
subject by: Kuhl et al. [26] where quasi-brittle materials where considered; Stumpf et al. [27] where the 
crack analysis based on the concept of continua with microstructure and evolving defects was discussed; 
Germain et al. [28] where the anisotropic layered material was formulated; Abu-Al-Rub et al. [29] were 
coupled anisotropic damage and plasticity constitutive model to predict the concrete distinct behavior in 
tension and compression was considered; Alastrue et al. [30] where fully three-dimensional anisotropic 
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elastic model for vascular tissue modelling was shown; and finally by Perzyna [31] and Sumelka et al. 
[32, 33] where class of implicitly non-local (rate type [34]) anisotropic models for metallic materials 
were considered. Nevertheless from the point of view of subject of this paper, the fractional anisotropic 
non-local models are nowadays still under development. 

In this paper we propose a generalisation of an original concept of isotropic fractional continuum me¬ 
chanics presented in [35, 36] for description of anisotropic non-locality. As will be discussed this for¬ 
mulation abandon not-only classical postulate of local action, but also restriction imposed by objectivity 
postulate (cf. discussion on objectivity in [37]). In this sense, this result will be analogous to the one 
obtained by Drapaca and Sivaloganathan [38], but is should be emphasised that both models can not be 
reduced one to each other, and operate in different physical dimension space. 

As mentioned the subject of this paper, namely the definition of non-local model utilising fractional cal¬ 
culus, belongs to the recent trends in mechanics [23, 24, 25, 39, 40, 41, 38]. In comparison to previous 
formulations the main advantage of the proposed formulation is its analogical structure to the general 
framework of classical continuum mechanics. Therefore, it allows, to some extent, to give similar phys¬ 
ical and geometrical meaning of introduced objects according to their classical counterparts (including 
of course some open questions due to interpretation of fractional calculus itself [42]). Other crucial 
advantages are: (i) we deal with finite deformations; (ii) the generalised fractional measures of the de¬ 
formation e.g. fractional deformation gradients or fractional strains have the same physical dimensions 
as classical one; (iii) characteristic length scale of the particular material is defined explicitly. 

The paper is structured as follows. In Section 2 existing fractional models of mechanics are summarised. 
The non-local fractional model accounting for anisotropic non-locality, is presented in Section 3. In 
Section 4, a benchmark examples are shown. 

Remark 1. Naming convention. Throughout the paper, we follow the naming convention common for 
the classical continuum mechanics. To show the influence of fractional calculus, we add the word frac¬ 
tional’. Therefore, we introduce objects such as fractional deformation gradient or fractional strain. 

It is important to notice, that similar nomenclature (because of fractional calculus application) 
is utilized in theories dealing with fractal media, where fractal/fractional continua is considered 
[43, 44, 45, 46, 47, 38, 48]. Those two concepts should not be confused. We read the influence of frac¬ 
tional differential operator, as in presented formulation, from phenomenological point of view leaving 
other interpretations for the future work. In other words, we examine some smoothed picture, without 
going into detail about the sub-scale constituents/phenomena. 


2 Existing fractional models of mechanics 

Let us shortly discuss selected concepts of non-local fractional continuum mechanics. 

Non-locality in the fractional continuum mechanics comes from the application of fractional derivative 
which operates over the interval. This interval defines simultaneously the range of non-local interaction. 
The fractional derivative means the derivative of an arbitrary order [49, 50], and in a special case when 
the applied order becomes integer we obtain classical local formulation (classical local derivative where 
the definition is given in a single point). In this sense, and because there are many different fractional 
differential operators [51], and many different ways how to introduce them into the specific formulation, 
we have nowadays several fractional continuum mechanics concepts. In general, one can distinguish 
two ways of introducing spatial non-local effects through fractional calculus into continuum mechanics 
[52]: (i) redefinition of kinematics; (ii) redefinition of constitutive law. 
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One of the first papers dealing with fractional kinematics was by Klimek [23] where the symmetric 
fractional derivative had been used in the definition of strain, namely 


1C a (x, t) = i ( 0 D° - x Dl) u(x, t ), 


( 1 ) 


where K. a is a fractional strain, D a is a fractional derivative, x is a spatial coordinate, t denotes time, 
and u is a displacement. In is important that the definition by Eq. (1) describes one dimensional problem 
under small strain assumption, the terminals 0 and L show that for calculation of fractional strain at the 
specific point of interest x we take all information from the body, and the fractional strain has physical 
. Another proposition was by Lazopoulos [25]. He proposed a measure of fractional strain 


1 —a 


unit [m 
in a form 


£ a {x, t ) = ( o-D" + x D a L ) t), 


( 2 ) 


where C a stands for fractional strain, with similar comments as for Eq. (1). Later a closely related 
proposition was stated by Atanackovic and Stankovic [40], namely 


£“(M) = \ ( -oc D x - x D Zo)u(x,t), 


(3) 


where £ a stands for fractional strain. 


Seminal work dealing with the proposition of fractional kinematics at finite strain for three dimensions 
was the paper by Drapaca and Sivaloganathan [38]. They started from the redefinition of motion, namely 


x = 


Kx aW x(X,t), 




—oo < an{t) <1, I,i = 1, 2, 3, 
an(t) = 1, I,i = 1,2,3, 


(4) 


where x is a motion, x is a spatial coordinate, X denotes a material coordinate, and a is an order of 
motion. In Eq. (4) 


K 


l-Q(t) 


K 


l-CCl(t) 

0 

0 


K 


0 


K 


0 

0 

i-as(i) 


(5) 


where ol { = (an, a 2 i, a 3i ), so 


a = 


a 11 Oti2 «13 
a 21 a 2 2 ct23 

Oi 31 «32 «33 


( 6 ) 


and finally 


Xi = 


_ _ pi — au pi— OL2i pi — _ 

— Xi — xi Xl u x 2 k x 3 x* — 


1 


r(l - aij(t))r(l - a 2 j(/))r(l - a 3 »(t)) 

_Xi(Y 'QdYidYidYa _ 

H \Xi- Yi !“«(*) |x 2 - y 2 \ a ^)\x 3 - y 3 | asi (0 : 


(V) 


where H = [ai, b\] x [a 2 , b 2 \ x [a 3 , 63 ] denotes the interval of non-local interaction, and T is the Euler 
gamma function 

roo 

T(a) = / e-H a - l dt. ( 8 ) 

Jo 
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Next, based on non-standard definition of motion by Eq. (4) the authors define the deformation gradient 
F a as 




dx(X,t) 

dX 


or 


pa _ 

b aA ~ 


d 


dX A 


( r \ 


—oil, 

X! 


R l x 3 °' 30 


Xa 


E, 


(9) 


where e a and E/i are base vectors in coordinate systems {x a \ and {X A }, respectively. 

Finally, the operator F a , which is now non-local, plays the same role as the classical deformation 
gradient. It is important to point out that F a has a physical unit, in this sense fractional strains also, 
and also for the first time the non-local fractional kinematics is defined on an interval which represents 
a part of the body, and enables to analyse anisotropic non-local effect. 


Concerning the second group of non-local fractional models, namely the one which bases on redef¬ 
inition of constitutive law, let us first mention about paper by Di Paola et al. [39]. In this paper the 
non-local elasticity relation was proposed as 


a(x, t ) — £ a £®cr(x, t ) = Es(x, t), 


( 10 ) 


where cr is a Cauchy stress, £ a is a length scale, £“ is a fractional derivative, and E denotes Young 
modulus. It is clear that Eq. (10) is a generalisation of well known Eringen model [53]. Later, Carpinteri 
et al. [41] considered a material constitutive law in a form 


a(x, t) = E 


f x + ^(C D> _c Di 


( 11 ) 


where k q is a material constant and has anomalous physical dimensions. Interestingly, Eq. (11) reflects 
the same basic formula as the one derived from Atanackovic and Stankovic [40] but starting from 
fractional non-local strain measure [54]. 

As a concluding remark, it should be emphasised that an interesting direction of research is also an 
attempt to unify non-local fractional models with theories dealing with fractal media or peridynamic 
model [55]. 


3 Fractional kinematics, stresses and balance laws 

3.1 Riesz-Caputo fractional derivative 


There are many definitions of fractional derivative [56, 49, 50, 57, 51]. They all share one common 
attribute i.e. they are all defined on an interval contrary to integer order differential operators defined 
in a point. In following part of this paper the Caputo’s type derivative over the interval (a, b ) is con¬ 
sidered. We call such operator Riesz-Caputo (RC) derivative cf. [58] and its definition includes linear 
combination of left and right Caputo’s derivatives, namely 

fwi = ^ ( Sivm+ (-D” ?»”/(<)). (i2) 

where ?(a) represents scalar valued function, is a left-sided Caputo’s derivative given by 

(t > a and n = [a] + 1) 


C 

a 


D?m 


i 

T(n — a) 



f (n \r) 

(■t - r) a ~ n+l 


dr, 


(13) 
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f D£f(t) is a right-sided Caputo derivative given by (t < b and n = [a] + 1) 




T(n - a) (t - t) a n+1 


(14) 


In the remaining part of this paper the RC derivative defined in Eq. (12) is shortly denoted as D a 
with the possibility of writing variable under the D in case of a partial differentiation of multivariate 

functions. For example D ° f represents partial fractional derivative of / with respect to the variable 

-X'l 

X\ over the interval which should be explicitly defined before X\ e (a, b). It is clear that for a = 1 we 
have 

™Dlf(t) = ^ f(t ). (15) 

Remark 2. Caputo’s derivative of a constant function. In general fractional derivative of a partic¬ 
ular type of a constant function is not equal zero. Caputo’s derivative is an exception, thus makes its 
application similar to classical integer order operator. 

Remark 3. On initial and boundary conditions. It is fundamental that using Caputo’s type derivative, 
one requires standard (like in the classical differential equations) initial and/or boundary conditions, 
whereas for other types of fractional derivatives (e.g. RL or GL) the initial/boundary conditions are of 
different type dependently on chosen definition. 

Remark 4. On the type of applied fractional differential operator. The type of applied fractional 
differential operator describes the type of non-locality. In other words, it governs the way in which 
the information from the surrounding influences particular point of interest. Thus, one should choose 
appropriate fractional differential operator, dependently on the material considered (in contrast the 
choice ofRC in the presented paper is due to reasons stated in Remarks 2 and 3). 

Recall, as an example, that the definition of the classical Riesz-Feller [59] fractional operator has an 
origin in processes with Levy stable probability distribution. In this sense it should be possible to define 
fractional differential operator in a way that e.g. it has an information about the distribution of grains 
sizes in a particular metal. 


3.2 Fractional deformation gradients 

It is important to emphasise that the specific fractional model presented in this paper was discussed in 
a series of papers devoted to fractional elasticity [60, 61]; fractional thermoelasticity [36], fractional 
Kirchhoff-Love Plates [62], and non-local rate independent plasticity [63]. Nevertheless, these results 
comprise the case of isotropic non-locality and it was shown that for such a case objectivity restriction 
known from classical mechanics holds. 

As mentioned in the introduction, in this paper we analyse the case of anisotropic non-locality i.e. 
materials showing scale effect but depending on direction. Therefore, this formulation abandon not- 
only classical postulate of local action, but also restriction imposed by objectivity postulate - what can 
be important for some classes of materials [37]. In this sense, this result is analogous to the one obtained 
by Drapaca and Sivaloganathan [38], but is should be emphasised that both models can not be reduced 
one to each other, and operate in different physical dimension space. 

The description is given in the Euclidean space. We refer to B as the reference configuration of the 
continuum body while S denotes its current configuration. Points in B are denoted by X and in S by x. 
Coordinate system for B is denoted by {Xa} with base E a and for S we have {x a } with base e a . 
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The regular motion of the material body B can be written as 

x = 0(X, t), 


(16) 


and its inverse as 

X = ¥>(x, t), (17) 

thus <f>t : B ->• S is a C 1 actual configuration of B in S, at time t. Here, it is important to emphasise that 
to avoid confusion we introduce new symbol for motion marked in Sec. 2 by x where we have followed 
Drapaca and Sivaloganathan [38] notation. In this sense, in the presented formulation the position of the 
non-local body is not changed contrary to [38], but the effort of the body (strains/stresses) is of course 
different due to non-local action. 

We define the fractional deformation gradient and its inverse as follows (cc e (0,1)) 

F(X, f) = £ a ~ 1 D° l (f)(X.,t), or F aA = i D a ^0 a e a ® E A , (18) 

XXX XX Xa 


and 


F ( X ,t) = l a ~ 1 D a cp( x ,t), or F Aa = £%*- 1 D a *><p A EA 


_ 0 OiAa — 1 


(19) 


where D a is a fractional differential operator in the sense of RC defined in previous section, and lx 

and l x are length scales in B and S, respectively. As an example the matrix representation of object F 

X 

is 


F = 

x 


’ l 

an¬ 

il 

_1 D ai1 

01 

l 

“12- 

12 

77“ 12 

01 

£ 

“13- 

13 

"t JJ&13 

01 

X 

X ! 

X 

X‘2 

X 

x 3 

l 

“21- 

21 

"I Y)OL2\ 

02 

l 

“22“ 

22 

-1 JJ&22 

02 

l 

“23" 

23 

"I ]JOi23 

02 

X 

Xx 

X 

x 2 

X 

X 3 

l 

“31- 

31 

-t 2)“31 

03 

l 

“32“ 

32 

"I Y) “32 

03 

l 

“33" 

33 

1 ]ja 33 

03 

X 

Xx 

X 

x 2 

X 

x 3 


( 20 ) 


Hence, in general we can write 


FF I = Sab^a ® E b, 

x X 


( 21 ) 


and 


FF i — fiab^a ® Gfe, 

X x 

where S tJ denotes the Kroncckcr delta. 

Having in mind classical definition of material dX and spatial dx line elements, namely 

dx = FdX, or dx a = F aA dX A e a , 


and inverse transformation 


dX = F x dx, or dX A = F\ 1 dx a E A , 


where 


dcj)CX.,t) d (j> a 

F(X, t) = ——, or F aA = T^e a ® E A 


F (x, t) = 


dX. 

dyX t) 

dx ' 


dX A 

or F Aa, = a^-F a ® e a . 


dx„ 


( 22 ) 

(23) 

(24) 

(25) 

(26) 
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we postulate the existence of fractional material and spatial line elements as 


dx = FdX, 

or 

dx a ■ 

— F' aA^X-A^ai 

X 

(27) 

and 





dX = Fdx, 

X 

or 

dX4 

— F Aa^-^aFiA’ 

X 

(28) 

Using Eqs (23),(24),(27) and (28) we have 





dx = FdX, 

or 

d x a 

cx. 

= F aA^XA^ci') 

(29) 

dX = FdX, 

X 

or 

dX B 

= F ba<1XaFb, 

X 

(30) 

dx = Fdx, 

X 

or 

d Xb 

a 

— F ba^Xa^b’) 

X 

(31) 


OL _ ^ Oi _ Q _ 

where F = FF _1 F _1 , F = FF and F = FF' 1 . Figure 1 summarises the situation. 

x x x x xx 


a 

F 



Figure 1: The relations between material and spatial line elements with their fractional counterparts 

Remark 5. On the non-locality of fractional deformation gradients. Classical deformation gradient F 

is local while fractional ones F and F are non-local due to the definition of Refractional differential op- 

X x 

erator. RC is based on an inten’al whose length depends on a particular material being described. Due 
to the definitions Eqs (18) and (19), their physical/geometrical interpretation, is analogous, compared 
with the classical one - the difference is that they operate on fractional or mixed (classical/factional) 
line elements. 
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Remark 6. On the length scale parameters in fractional deformation gradients definitions. The main 
attention should be paid to the length scale parameters appearing in definitions Eqs (18) and (19) of 
the fractional deformation gradients. 

It is important to notice that without those parameters, the unit of the fractional deformation gradients 
would be (in SI) rn l ~ a . Therefore, the introduction of the length scales, similarly to classical gradi¬ 
ent continuum models, allows to finally obtain dimensionless quantity. Therefore, we can compare the 
lengths of line elements dX and dx with fractional ones dX and dx what is crucial concerning possible 
strains definitions. 

In previous papers, where isotropic non-locality was considered as mentioned, it was shown also, that 
it is necessary to introduce the length scales in order to fidfil the rigid body motion requirements (ob¬ 
jectivity restriction). Nevertheless, from purely mathematical point of view, those parameters could be 
omitted. 


Remark 7. On the equivalence of the fractional formulation with the classical one for <x %3 = 1. For 

= 1 RC continua becomes classical one and hence we have: 


a. 


v 


F = F = F 1 = F. 

X X 


(32) 


F^ 1 = F _1 = F = F 

x x ’ 


a 

F — I, 




dx = dx, 


(33) 

(34) 

(35) 

(36) 


dX = dX. 


(37) 


Remark 8. On the fractional deformation gradient properties. F and F are not two point tensors 

x X 

while F, F and F are. Based on the properties of motion Eq. (16) we have that the inverse of F and F 

x x x x 

exists. 


3.3 Fractional strains 

We define the strains by analogy to the classical continuum mechanics based on the difference in scalar 
products in actual and reference configurations. The introduced fractional deformation gradients allows 
to define 4 concepts of strains. Thus, one can define: 


E=-(F 7 F — I), or Eab = A a F aB — Iab)Ea®E,b, 


(38) 
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1 o T <> 1 o o 

e= 2^ i_F F ^ ° r 6ab = 2( iab ~ F «A F ^ ea ® eb ' 


(39) 


where E is the classical Green-Lagrange strain tensor or its fractional counterpart; e is the classical 

o 

Euler-Almansi strain tensor or its fractional counterpart, and depending on the formulation, F can be 

OL O 

replaced with F or F or F or F. According to the chosen F the associated others variables like (they 

X x 

all have, to some extent, classical meaning according to Remark 5, but simultaneously are non-local): 

the left and right Cauchy-Green tensors; the orthogonal tensor and left or right stretch tensor from polar 

o 

decomposition of F can be defined using analogical to classical rules. 

Remark 9. On the relation between fractional displacement gradient tensor and fractional strains. 

As in the classical continuum mechanics one can define the relation between strains and displacement 

gradient tensor utilising introduced fractional gradient tensors F and F. 

x x 

The displacements in the material description U are defined as: 


U(X,i) = x(X, t) — X, 


(40) 


and its fractional gradient 

GradU = F — I, or l “^" 1 D a ^U a = (F aA 

xx x a/1 x A x 


thus we have 

F = GradU + I 

x x 

Similarly, the displacements in spatial description u are defined as 

u(x, t) = x — X(x, t), 


hi /I ) C„ 0 E^4, 


(41) 

(42) 


(43) 


and its fractional gradient 

gradu = i - F, or i ^ a_1 D aAa u A = (Ua - F Aa) Ea 0 e a , (44) 

thus we have 

F = i — gradu. (45) 

X X 


By applying Eqs (42) and (45) into the fractional strain definitions Eqs (38) and (39) we obtain their 
dependence on the fractional displacement gradients. 


As in classical continuum mechanics, we can introduce small fractional Cauchy strain tensor, we have 

( 1 = 1 = 1 ) 

X x 


o 1 

£ = - 
2 


gradu + gradu 2 


o 


where £ stands for fractional Cauchy strain. For ID deformation Eq. (46) reads 


| i ( x _aD%u - C x D x+b u) . 


(46) 


(47) 


Notice, that Eq. (47) is similar to definitions presented in [23] (Eq. 9), [25], [40] (Eq. 2.7) and [41] 
(Eq. 20), but operates on finite internal (similarly to the ’short memory’ principle discussed in [49]), 
and gives non-dimensional quantity. It should be emphasised, that in contrast to isotropic non-locality 
discussed in previous papers (cf. [36]) a f b, but it remains physically reasonable to assume some 
relation between £ and length of internal over which fractional derivative is calculates e.g. equivalence. 
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3.4 Stresses and balance laws 


According to Fig. 1, material line element is dX or dX, while spatial dx or dx (together with cor¬ 
responding fractional/classical deformation gradient). The selection of a specific path causes that the 
remaining one can be considered utilising the concept of dual variables [64], For example, assuming 
that we describe the deformation in terms of dX and dx, classical deformation (based on dX and dx) 
appears as (one of possible) ’intermediate’ one. Therefore, it is clear that stresses, and balance laws 
should be postulated dependently on chosen description of deformation in a standard manner (such 
result is analogous to the one presented in [40] Eqs 2.7,2.17, and 2.18). 

To understand this logic (without loss of generality) let us consider purely mechanical problem. It will 
be observed that the restrictions to be held by fractional stresses are analogous to the classical one, but 
fulfilled in auxiliary ’fractional’ (or phenomenological) space. 

First, following Fig. 1, the material (spatial) volume element dF (dn) or its fractional counterpart 
dV (dn), and the material (spatial) surface element dS (ds) or its fractional counterpart dS (ds) are 
introduced - cf. [60]. Next, the fractional Cauchy (true) traction vector t(„) exerted on ds with outward 
normal n is obtained as a transformation of classical Cauchy (true) traction vector t( n ) exerted on ds 
with outward normal n. Thus, according to the diagram in Fig. 1 we have 

_ a. a. 

t(n) Or i(n) b F bat(n) a^b- (48) 

As in classical set-up we postulate the relationship between fractional traction and fractional Cauchy 
stress tensor 

t(n) = nor, (49) 

where a denotes fractional Cauchy stress tensor. Based on the above relations, and the fact that 

n = Fn, (50) 

X 


we have 

a a ex. 

Ft/' n 'i = Fner = Fncr, 

X y ’ X X 

thus finally 

cr = a, 


(51) 

(52) 


where a is classical Cauchy stress tensor. Thus, according to relation Eq. (52) the fractional stresses 
are transformation of the classical Cauchy measure to auxiliary ’fractional’ space. 

Finally, before the balance of momentum in spatial description is postulated for fractional case, the 
conservation of mass is considered in a form 



p$dV = podV = pdv = pdv, 

(53) 

or shortly 

o ^.o 

Po = J Pi 

(54) 

where po (po) is the reference mass density (fractional counterpart), p (p) is 

o o 

tional counterpart), and J = detF is a Jacobian (as before in the denotation 
classical quantity or fractional counterpart). The following relation holds 

spatial mass density (frac- 

0 

(•) one can replace (•) with 


f o? .0 /■ <>.0 f 0 0 

Ipvdv = ltd s+ Ipfdv, 

J v J dv J v 

(55) 
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where v is a velocity, and f is a body force per unit mass. By applying the divergence theorem to 
Eq. (55) we have 



or in the absence of inertia forces 

O t O 

diver T + pf = 0. (58) 

Furthermore, the symmetry of a (so er or <x) can be checked in classical manner utilising the balance 
of moment of momentum. 


As a concluding remark, it is clear that formulation in a material description needs proper definitions 
for the first and the second (fractional) Piola-Kirchhoff stress tensors. Thus, dependently on chosen 
fractional strains definition we have: 


p 

= JcrF~ T , 

S 

= 

JF-'crF- 7, = 



~ -T 




p 

= JcrF . 

s 

= 

JF crF = 

X 


X 

X 


XX X 



1 ~ -T 



~_1~ ~-T 

p 

= r 

a F 

i 

s 

= J FcrF 

X 

X 

X 


X 

XXX 

a 

a 

a-T 

a 


a a ~ 1 a 

p 

= JcrF , 

s 

= 

,JF crF = 


F -1 P, (59) 

F _1 P, (60) 

x x 

= FP, (61) 

X X 

a — 1 a 

F P, (62) 


_ _ a _ _ a 

where P, P, P, P and S, S, S, S are classical/fractional the first and the second Piola-Kirchhoff stress 

x x x x 

tensors, respectively. As previously, taking = 1 (local kinematics) the introduced stresses definitions 
reduce to the corresponding classical one 


P = P = P = P, (63) 

X x 

s = s = S = S. (64) 

X x 

Remark 10. On physical interpretation of the parameters a and i. Introducing non-locality utilising 
fractional calculus we simultaneously add a set material parameters a and £. The order of fractional 
continua a {together with type of applied fractional differential operator) controls the way in which the 
information from the surrounding influences particular point of interest. The length scale t defines the 
amount of this information {size of this surrounding). So, both parameters are crucial. In the following 
example we will observe that both a -> 1 and £ -> 0 recovers classical solution - as it should be. 
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4 Examples 

In the example section, the influence of fractional kinematics on the fractional strains is presented based 
on the assumed different motions. The following notations for strains are proposed (cf. Eqs (38) and 
(39)): 

1. Classical formulation 

E = 1(F t F-I), e = i(i —F^F- 1 ). (65) 

2. Formulation based on the fractional spatial line element (dx) and the classical material line ele¬ 
ment (dX) 


E = -(F F — I) 


x 


2 v x x 


e = -(i-F 
x 2 v x 


-T 


F - 1 ). 

X 


( 66 ) 


3. Formulation based on the classical spatial line element (dx) and the fractional material line ele¬ 
ment (dX) 


E = -(F -t F -1 

X 2 X X 




F t F). 


(67) 


4. Formulation based on the fractional spatial line element (dx) and the fractional material line 
element (dX) 





( 68 ) 


In following examples ID deformation is considered. We assume for brevity that: (i) l = l = l\ 

x x 

(ii) 0 < a < 1; (iii) the deformation in a particular point of interest X,\ is calculated for terminals 
a = Xa — f-L and b = £r + Xa (this is because of anisotropy of non-locality that Xa in general does 
not necessarily lay in the middle of interval (a, 6)). 


4.1 Example 1 

Fet us consider a particular type of motion, namely the case when motion is a linear function of co¬ 
ordinate. We compare for brevity the deformation based on F and F. Thus, the motion is described 
as 

x = </>(X) = (1 + /3)Xiei + X- 2 e 2 + X 3 e 3 . (69) 


For such a case we have 


F = 


(1 + /3) 0 O' 

0 1 0 

0 0 1 


and (we omit the underset X for <,) 


(1 + /3) 0 O' 

0 1 0 

0 0 1 


(70) 


F = M 

x 
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thus 


and 


E = - 
2 


'(1 + /3) 2 - 1 0 O' 

0 0 0 

0 0 0 


E = - 
x 2 


[A4(l + /?)] 2 -l 
0 
0 


0 

M 2 - 1 


0 

0 

M 2 - 1 


(72) 


(73) 


where M. = 1 (Z\~ a + 7^7 °j ■ We observe, that anisotropic fractional deformation intro¬ 

duces in general directional deformation e.g. due to evolving micro-structure or existence of electro¬ 
magnetic, thermic, and/or chemical processes [38]. 


4.2 Example 2 

Lets consider a motion being a non-linear function of the material coordinates Xa- As before, for 
brevity, the deformation based on F and F is compared. Thus, the motion is described by 


x = (f>(X) = e Xl ei + X 2 e 2 + X 3 e 3 . 
In this case the classical solution is presented as 



e Xl 

0 

o' 


~e 2X] - 1 

0 

O' 

F = 

0 

1 

0 

0 

0 

0 


_ 0 

0 

1 


0 

0 

0 


(74) 


(75) 


The calculation of F requires detailed explanation. By analogy to the first example F 22 = F 33 = M. 

X X X 

whereas F n will be approximated numerically, because analytical solutions under fractional calculus 

X 

are very limited [65, 57]. Therefore, according to Eq. (12) we need to calculate adequate left and right 
Caputo derivatives for (f>. We have 


F 

x 


11 — 


pot —1 RC 

' a 


D 

x 1 


_ 

b e ~ 



c n a p Xi 

X 1 -i L JJ x 1 e 


c 

X 1 


D 


a 

X x +t R 


a Xx 


)■ 


(76) 


where for calculations puiposes ?(a) = T (2 - a), and l = \{l L + t R ). 

The following approximations utilising the modified trapezoidal rule can be used for the calculations 
[57]. For the left sided derivatives we use: 

a = to < t\ < ... < tj < ... < t m = t, h = — - — = - m > 2, (77) 

mm 


un—a , 

f = r(n _ a + 2) {[(™ - - (m - n + a - + 

m— 1 

f {n) (t m ) + [(m - j + l) n - Q+1 - 2(m - j) n ~ a+1 +(m-j- l) n ~ a+1 ]f^(t,)} , (78) 

3 =1 

where f denotes classical n-th derivative at t = tj. 


13 
















Wojciech Sumelka - Compiled on 2015/02/09 at 01:27:28 


Similarly, for the right sided derivatives we use: 

t = to < ti < ... < tj < ... < t m = b, h = —-— = ^ m > 2 , (79) 

m m 

un—a . 

- (-1)"--—- { [(m - 1)""“ +1 -{m-n + a- l)m n ~ a }f^(t m ) + 

i (n — a + 2) t 

m— 1 

/ W (*o) + E [O' + - 2 ^' n ' a+1 + 0 - l)"-«+ 1 ]/W(t i )} • (80) 

7=1 

The obtained results are shown in Figs (2), (3) and (4), where the comparison of fractional strain mea¬ 
sure E against the anisotropy of non-locality, the order of derivative a, and length scale i is presented. 

X 

It should be noticed that for l —> 0 (thus, simultaneously the length of the interval over which fractional 
measures are calculated approaches to zero) fractional strain measure merge with the classical definition 
(which are defined in a point) independently of the chosen order of fractional derivative and anisotropy 
of non-locality. Similarly, if a —> 1 once more the classical strains arc captured, however independently 
of chosen l. 

Another interesting result is the observation that one could model size effect utilising fractional strains. 
It is clearly presented in Figs (2), (3) and (4) where the dependence of the results against assumed length 
scale t is shown. 

Recall, that the introduced fractional strains measures have meaning analogous to the classical one. 
The difference is that the particular value of strains in a point includes the information from its finite 
suiTounding (controlled by II and (/,>) in the way defined by order of fractional continua (controlled by 

a). 

5 Conclusions 

In this paper the new concept of generalisation of the classical continuum mechanics utilising fractional 
calculus is presented. Due to fractional derivative properties the obtained formulation is non-local, and 
this non-locality is of anisotropic type. It should be pointed out that the defined non-local model can be 
applied to any known constitutive model of elasticity, hyper/hypo-elasticity, visco/plasticity etc. 

The anisotropic non-locality introduces a set material parameters: orders of fractional continua (aij) 
and length scales (lij). They control the way in which the information from the surrounding influences 
particular point of interest. This result by far enhances modelling approach and allows closely mimic 
the experimental observations. 

Based on two examples picturing one dimensional deformation, the role of the non-local definition of 
strains based on the fractional differential operators is presented along with the numerical treatment of 
the fractional derivatives. Concluding, fractional models of mechanics give the opportunity for deeper 
insight into the analysed problem contrary to classical local formulation. 
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